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ABSTRACT 

By means of galaxy evolutionary models, we explore the direct consequences of the A Cold 
Dark Matter (ACDM) cosmogony on the size evolution of galactic discs, avoiding intention- 
I ally the introduction of intermediate (uncertain) astrophysical processes. Based on the shape 

of the rotation curves and guided by a simplicity criterion, we adopt an average galaxy mass 
baryon fraction /gai = 0.03. In order to study general behaviors, only models with the av- 
erage initial conditions are analyzed. The stellar and B-band effective radii, i?* and Rb, of 
O ' individual galaxies grow significantly with time (inside-out disc formation) with laws that are 

weakly dependent on stellar mass, M^, or luminosity. Lb- However, the change of i?, with z 
^ . at a^xeii Af, is slow; for z ^ 2.5, i?,(M*=constant)(x (l + z)"""' for a large range of masses. 

^ I ' On the other hand, the change of Rb with z at Si fixed Lb is strong and it resembles the R b de- 

creasing law of the individual models; roughly Rb{L B=constant)oc (1 + z)^'^-^^ for z < 0.75, 
CN ■ and oc (1 + z)''^-^ for higher z's. We find also that at z = 0, i?, cx Af,"'^^ and Rb oc L°g'^°, 

\ remaining the slopes of these relations practically the same up to z w 3. Our model predic- 

tions are in reasonable agreement with observational inferences on the typical radius change 
with z of late-type galaxies more luminous (massive) than high values imposed by the se- 
lection effects. The models seem also to be consistent, within the large scatter, with the Rb 
and Lb values obtained from small (non complete) samples of sub-L* late-type galaxies with 
I available rest-frame photometric information at different z's. The properties and evolution of 

the ACDM haloes seem to be the main drivers of galaxy disk size evolution. Nevertheless, the 
models reveal a potential difficulty in explaining the observed steepening of the Rb-Lb rela- 
tion with respect to the R^-M^, one, an effect related to the well established color-magnitude 
relation. 
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1 INTRODUCTION 

The completion of multi-wavelength deep field surveys during the 
last years has enormously benefited the study of galaxy evolution. 
The analysis of galaxies with redshift known from these surveys 
provides valuable information about the changes with cosmic time 
of galaxy population properties. However, since the evolution of 
individual galaxies is not directly observable and because of the 
strong selection effects at higher redshifts, a direct interpretation 
of the observations is not an easy task. An adequate comparison of 
theoretical models with observations helps largely in this undertak- 
ing. Here we will present evolutionary models of disc galaxy for- 
mation and evolution inside growing A Cold Dark Matter (ACDM) 
haloes with the aim to show predictions (i) on size evolution of in- 
dividual galaxies, and (ii) on the typical size change with redshift 
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z of galaxies of a fixed stellar mass, M,, or luminosity, L. The 
latter results are those that can be compared with the information 
provided by observations. 

The qualitative description of disc galaxy formation in ACDM 
haloes is based on the idea that both dark and baryonic matter ac- 
quire angular momentum by tidal torques while the perturbation 
is in its linear regime. The gas within the virialized halo cools 
and falls to its center until it attains centrifugal equilibrium form- 
ing a disc (White & Rees 1978). According to the spherical col- 
lapse model, an overdense region of mass A/h virializes with a 
radius _Rh at the epoch z when its cumulative density overcomes 
Pbg(z) X Ac{z), where pbs{z) is the background density at z and 
Ac is a critical value at each epoch z that depends on the cosmo- 
logical model (Bryan & Norman 1998). It follows then that 
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where H{z) is the Hubble parameter. If the disc galaxy radius 
is proportional to i?h, then a significant decrease of sizes with z 
at a fixed halo mass is expected (e.g., Mo, Mao, & White 1998; 
Bouwens & Silk 2002). For haloes treated as truncated isothermal 
spheres and assuming that an exponential disc is formed under de- 
tailed angular momentum conservation, it indeed follows that h^ cx 
Ai?h (e.g.. Fall & Efstathiou 1980; Mo et al. 1998), where hd is the 
disc scale length and A is the halo spin parameter (Peebles 1969). 
Taking into account a Navarro-Frenk- White (NFW) density profile 
and including the adiabatic halo contraction due to disc formation, 
the latter dependence becomes hd oc Ai?h }c ^^'^ c, /gai), 

where c is the concentration parameter of the NFW profile, fa is a 
shape-function that increases with c, /gai = A/bar /A#h is the mass 
galaxy baryon fraction, and fa is a function that takes into account 
the NFW halo adiabatic contraction for an exponential baryonic 
disc (Mo et al. 1998). Somerville et al. (2008) have shown that a 
further reduction of the change of galaxy radius with z at a fixed 
mass is obtained considering the fact that c decreases roughly as 
( 1 + «) ~ ^ . These authors concluded that the observed weak change 
with z of typical stellar radii of disc galaxies of similar stellar 
masses at all redshifts agrees with the model predictions under the 
key assumption that A and /gai do not change with z. 

The models based on the Mo ct al. (1998) approach (e.g., those 
of Somerville et al. 2008) actually refer to a 'static' (instantaneous) 
population of disc galaxies at a given z; the models do not follow 
actually the individual and local evolution of the halo-galaxy sys- 
tem. Besides, such models consider only baryonic discs and assume 
they have an exponential mass surface density distribution. The 
complex evolutionary processes of gas infall and local gas transfor- 
mation into stars, with the consequent prediction of gas, stellar, and 
luminous disc properties, are beyond the scope of the mentioned 
models. The cosmological N-body/hydrodynamical simulations are 
on the other extreme of complexity. These simulations use to report 
the so-called 'angular momentum catastrophe' with the consequent 
formation of too small and compact discs (Navarro & White 1994). 
Such a problem seems to have been largely ameliorated in recent 
simulations of small volumes with a very high resolution and an ap- 
propriate tuning of sub-grid physics (e.g., Abadi et al. 2003; Gov- 
emato et al. 2004,2007; Robertson et al. 2004; Zavala, Okamoto & 
Frenk 2008; Scannapieco et al. 2008). Even so, the full numerical 
simulations do not yet allow to predict in detail the evolution of a 
whole population of well resolved galaxies. 

From the observational side, the study of galaxy size evolu- 
tion is hampered by the flux and surface brightness (SB) limits, 
which introduce incompleteness and selection effects in the pho- 
tometric samples, particularly for those at higher redshifts. Thus, 
the observational inferences are yet controversial (see for recent re- 
views Cameron & Driver 2007; Elmegreen et al. 2007). Because 
magnitude and size correlate strongly, the determination of the av- 
erage rest-frame surface brightnesses (SB) as a function of z has 
been used commonly to parametrize the evolution of galaxy sizes. 
In this analysis is crucial the choice of the selection and complete- 
ness functions to be applied to each z bin. In early works, an in- 
crease of the mean rest-frame B— band SB of ~ 1 mag from z ~ Q 
to z ~ 1 was reported for late-type galaxies (Schade et al. 1995; 
Roche et al. 1998). Later on, the use of stricter completeness func- 
tions, defined for the highest z bins and imposed to the lower-« 
bins, showed no detectable mean SB evolution (Simard et al. 1999; 
Ravindranath et al. 2004). More recently, larger galaxy samples 
were used, and the completeness functions were estimated indepen- 
dently for each z bin by means of artificial galaxy simulations. As a 
result, a non-negUgible SB or radius (at a given luminosity) evolu- 



tion for luminous late-type galaxies was inferred (e.g., Bouwens & 
Silk 2002; Bouwens et al. 2004; Barden et al. 2005; Trujillo et al. 
2006; Cameron & Driver 2007; Franx et al. 2008). Some of these 
samples extend up to z ~ 3 or even z ^ Q. However, as higher is 
z, (i) only the most luminous and highest SB galaxies of the epoch 
are observed, and (ii) the determination of the galaxy type becomes 
very uncertain. 

It is important to have in mind that, likely, both radius and lu- 
minosity evolve. Therefore, is not straightforward to decipher from 
the sample mean SB change with z how is the size evolution of in- 
dividual discs. Besides, the SB evolution inferences happen to be 
different in different rest-frame bands. Additional kinematic infor- 
mation (rotation velocities) certainly may help on this undertaken 
(see e.g., Mao, Mo & White 1998). 

Given the current status of the theory and observations men- 
tioned above, the following natural questions arise. Do the direct 
implications of the ACDM cosmogony on galactic disc evolution 
agree with the observational inferences? Could galaxy evolutionary 
models shed light on the connection between the individual evolu- 
tion of discs and the sample properties inferred from photometric 
observations? 

With the aim to answer these questions, we will use here the 
semi-numerical evolutionary approach: the assembling of haloes 
and the discs within them is (i) followed self-consistently (for the 
discs we use hydrodynamic techniques and take into accoimt real- 
istic descriptions for the relevant astrophysical processes, e.g., star 
formation, turbulence and its dissipation, stellar dynamics, gravi- 
tational instability, etc.), and (ii) for a large range of initial con- 
ditions (Avila-Reese et al. 1998; Firmani & Avila-Reese 2000 - 
hereafter FA2000; Avila-Reese & Fkmani 2000,2001). Alternative 
approaches with a similar philosophy have been developed by other 
authors (e.g., van den Bosch 2000; Naab & Ostriker 2006; Stringer 
& Benson 2007; Dutton & van den Bosch 2008). As a result of all 
these approaches, the models show that discs evolve inside out, and 
the predicted baryonic, stellar and luminous present-day scaling re- 
lations tend to be in agreement with observations (see e.g., FAOO; 
Avila-Reese et al. 2008; Dutton & van den Bosch 2008). Here we 
will explore the direct predictions of such models regarding disk 
size evolution. With this aim in mind, we will not include the possi- 
ble effects of intermediate (uncertain) astrophysical processes [e.g., 
mass outflows, long time gas cooling, active galactic nuclei (AGN) 
feedback]. 

In Section 2 a brief description of the galaxy evolutionary 
models and the strategy followed in this paper are presented. The 
results regarding disc size evolution and comparison with observa- 
tional inferences are given in Section 3. Finally, a summary and the 
conclusions of the paper are presented in Section 4. We assume a 
concordance cosmology with h = 0.7, Oa = 0.7, flm = 0.3, 
«(, = 0.04, and as = 0.8. 



2 THE MODEL 

Our models calculate the formation and evolution of disc galax- 
ies inside growing ACDM haloes, attempting to keep the underly- 
ing physics as transparent as possible. The main physical ingredi- 
ents are as follows (see for details FA2000; Avila-Reese & Firmani 
2000). An extended Press-Schechter approach is used to gener- 
ate the halo mass accretion histories (MAHs) from the primordial 
Gaussian density fluctuation field. A generalized secondary infall 
model with elUptical orbits is applied to calculate the time-by-time 
viriaUzation of the accreting mass shells (Avila-Reese et al. 1998). 
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Figure 1. Total rotation curves for "central" models of Afjj = 7 10^ Mq 
(bottom green) and 7 IO^^Mq (top red). From top to bottom, the curves 
correspond to z =0, 1, and 3. Regions smaller than twice the effective 
radius (approximately one optical radius) are shown with the thick trace. 



The orbit ellipticity is fixed in such a way that the structure of the 
ACDM haloes agrees well with results from cosmological N~body 
simulations (Avila-Reese et al. 1999; FA2000). 

A constant (baryon) fraction /gai of the mass of each shell is 
assumed to cool down in a dynamical time and form a disc layer. 
The accreting mass shells at the time of their virialization, t„ , are 
assumed to rotate rigidly with a specific angular momentum calcu- 
lated as jsh{tv) = AJh/AMh, where A represents a difference 
between two time steps, and A = XhGM^^'^ / \Eh\^^'^, Mh, and 
Eh, are the halo total angular momentum, mass, and energy re- 
spectively; Aft is the halo spin parameter, assumed to be constant 
in time. As a result of the assembling of these mass shells, a present- 
day halo ends with an angular momentum distribution close to the 
(universal) distribution measured in N-body simulations (Bullock 
et al. 2001). The radial mass distribution of the layer is calculated 
by equating its specific angular momentum to that of its final cir- 
cular orbit in centrifugal equilibrium (detailed angular momentum 
conservation). The superposition of these layers form a gaseous 
disc, which tends to be steeper in the centre and flatter at the pe- 
riphery than the exponential law. The gravitational interaction of 
disc and inner halo during their assembly is calculated using the 
adiabatic invariance formalism (Blumenthal et al. 1986). 

A further step is the calculation of the stellar surface density 
profile. The previously mentioned processes and the fact that star 
formation (SF) is less efficient at the periphery than in the centre, 
produces stellar discs with a nearly exponential surface density dis- 
tribution, and the size mainly determined by A^. A cusp is present 
in the inner region, where the bulge is identified in observed galax- 
ies. The disc SF at a given radius (assuming azimuthal symmetry) is 
triggered by the Toomre gas gravitational instability criterion and 
self-regulated by a balance between the energy input due to SNe 
and the turbulent energy dissipation in the ISM. Both such ingre- 
dients determine the gas disc height and the SF rate. This physical 
prescription naturally yields a Schmidt-Kennicutt-like law. The SF 
efficiency depends on the gas surface density determined mainly 
by \h, and on the gas infall history, which is assumed to be propor- 
tional to the halo MAH. Very simple population synthesis models 
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Figure 2. Dotted curves are the _R» (upper panel) and Rg (lower panel) 
evolutionary paths of galaxy models that at z = have halo virial masses 
(from top to bottom) Mh/MQ= 2.8 10^^, 7.0 lO^^, 4.4 lO^^, 2.8 lO^^, 
1.76 10^2^ I I iQi2^ 7 Q ]^Qii_ 4 4 ;^Qii_ 2.8 10i\ 1.76 IQ^'^ , 1.1 lO", 
and 7.0 10^" (logarithmic steps of 0.2). The solid curves connect points of 
constant LogA/» and LogLg, respectively along the individual models. The 
corresponding constant values are indicated inside the panels; the unities 
are Mq for M« (upper panel), and Lbq for Lb (lower panel). Symbols 
in the lower panel correspond to observed galaxies (Tamm & Tenjes 2006): 
hogLg <9.5 (red triangles), 9.5 ^LogL^ ^ 10.0 (green squares), and 
LogLg > 10.0 (blue circles). 



are used to calculate B— and band luminosities from the SF 
rates. 

A given galaxy model is defined then by the halo MAH and 
Xh, and by /gai- The idea is to keep the model as simple as possi- 
ble in order to explore the direct consequences of the ACDM cos- 
mogony on disc galaxy evolution. Since we are interested here in 
the generic evolutionary trends, we study only the "central" mod- 
els of different masses characterized by: 

• the averaged MAH corresponding to the given halo mass Afh, 

• a value of Ah — 0.035, which is close to the mean of relaxed 
haloes measured in numerical simulations (e.g.. Belt et al. 2007), 

• an average value of /gai = 0.03 

Note that all these physical ingredients have in fact dispersions (not 
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taken into account here) that, of course, will produce a wide scatter 
in the values of the output galaxy properties. 

A remark has to be made about /gai. i^everal pieces of evi- 
dence suggest values of /gai < 0.05, though /gai is expected to 
vary with mass, being smaller for both very low and very high lumi- 
nosity galaxies. Some physical mechanisms suggested for produc- 
ing such low values are SN-driven galaxy mass outflows, long gas 
cooling times in massive haloes, SN and AGN halo gas re-heating. 
Given the current uncertainties in the understanding of these pro- 
cesses as well as in the observational determination of the halo 
masses, Mh, we have preferred here to use a constant value for 
/gal equal to 0.03. With such a value, our disc galaxy models agree 
with observed dynamical properties, mainly the approximate flat- 
ness of the rotation curves at different epochs and the present-day 
disc-to-total velocity ratios at the maximum of the rotation curve 
(see e.g.. Mo et al. 1998; FA2000; Zavala et al. 2003; Gnedin et al. 
2007; Dutton et al. 2007). 

In Fig. [T] we plot the total rotation curve evolution {z — 
0,1 and 3) corresponding to our "central" I models with Mh = 
7 1O^°M0 (bottom green) and 7 lO^^M© (top red). The curves are 
nearly flat around twice the effective stellar radius _R, (indicated 
by the thick trace), where R, is the radius where half of the total 
stellar mass Af* is contained; two R, correspond approximately to 
one optical radius. 



3 RESULTS 

In Fig.m the dotted lines show the evolution of the stellar effective 
radius, i?, (upper panel), and the B— band effective radius, Rb 
(lower panel) for models with total present-day virial masses rang- 
ing from A4 = 7 IO^'^Mq to A4 = 2.8 IO^^Mq by steps of 0.2 in 
logarithm. Both 7?* and Rb significantly decrease with z, the de- 
creasing being slightly steeper for the more massive galaxies. For 
present-day disc stellar masses ranging from ~ 10® to ~ lO^^M©, 
the stellar effective radii decreases up to 2 = 1 by factors from 1.75 
to 1.9, respectively; up to 2 = 2, these factors are from 2.8 to 3.2, 
respectively. In the B band, for the same models, which at z = 
have Lb from ~ 4.2 10^" to ~ 5.3 10* Lbq, the corresponding 
radius decreasing factors up to 2 = 1 are from 1.75 to 1.85, and up 
to 2 = 2, are from 2.8 to 3.3, respectively. The stellar and B-band 
radius evolutions tend to be similar on average, while Rb tends to 
be 1.4-1.6 times larger than R, (at 2 — 0). 

Along each curve in Fig. (2] Al, and Lb are also changing. 
We may connect the points in the different curves that have the 
same M, value (upper panel) and Lb value (lower panel). These 
curves (solid line) show the radius that disc galaxies of a fixed M* 
(Lb ) have at different epochs; this is actually the kind of informa- 
tion that obsen'ations can provide, rather than the individual ra- 
dius evolution. The change of radius with z along the A/*=constant 
curves is small. For z < 2.5, i?,(A/*=constant) decreases roughly 
as (1 + z)"" * for all the masses studied here. In the case of the 
L_B=constant curves, the effective radius decrease with 2 is strong 
and it resembles the decreasing law of the individual models. For 
z ^ 0.75, 7?s(L_B=constant) decreases roug hly as (l + 2)-''■*^and 
for 0.75 S z < 3 it does as (1 + 2)"^ ^ 

The different radius dependences on 2 for models in the upper 
and lower panels of Fig.|2]are explained basically by the differences 
in the evolution of M* and Lb, which implies a significant change 
with 2 of the M^ILb ratio. These different evolutionary paths can 
also be appreciated in the R,{z) vs Mt{z) and Rb{z) vs Lb{z) 
diagrams of Fig. [3] (upper and lower panels), respectively, where 




Figure 3. Relations R^-Mt (upper panel) and Rb-Lb (lower panel) 
at four epochs, indicated inside the panels as Log(l+2), for the "cen- 
tral" model galaxies studied here. Dashed lines are the evolutionary tracks 
of individual models. Symbols in the lower panel correspond to ob- 
served galaxies (Tamm & Tenjes 2006): Log(l+2)< 0.3 (blue triangles), 
0.3 ^Log(l+z)5C 0.5 (green squares), and Log(l+2:)> 0.5 (red circles) 



the solid lines connect models at a given 2, while the dashed lines 
show the evolutionary tracks of individual models. At least for the 
"central" models studied here, we see that at z = 0, i?* oc MJ^'"^^ 
and Rb oc L'^b'^^\ Furthermore, the slopes of these correlations 
roughly remain the same at higher redshifts in such a way that the 
changes with 2 of the relation zero-points, i?, (2)/!/, (2)"'^** and 
Rb{z) / LB^zf''^'^ , describe basically the curves of radius change 
at M»=constant and Ls=constant. plotted in Fig.[2] Thus, our mod- 
els show that using the radius change with 2 for a fixed mass (lu- 
minosity), or using the zero-point evolution of the radius-mass (- 
luminosity) relations for inferring the size evolution of galaxies, 
depends on how A/, and L_g do evolve. 

Our results agree conceptually with those by Somerville et al. 
(2008). However, the predictions of these authors refer to 'static' 
populations and to baryonic discs rather than to individual evolving 
stellar and luminous discs (their approach does not allow to calcu- 
late gas transformation into stars). With our evolutionary models, 
we have generated a plot similar to the upper panel of Fig.|2]but for 
the baryonic effective radius and for curves of Afbar=constant. The 
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obtained curves follow very closely the stellar ones, showing that 
the radius evolution of the baryonic and stellar disc is in fact quite 
similar. 



3.1 Model results vs observational inferences 

In the lower panels of Figs. |2] and [3] the data corresponding to 22 
disc galaxies with photometric z's between 0.5 and 2.6 reported in 
Tamm & Tenjes (2006) are plotted. The authors used deep Near 
Infrared Camera and Multi-Object Spectrometer (NICMOS) J and 
H band images from the HDF-S field to infer rest-frame B-band 
photometric properties. The sample is not complete in any sense 
but it gives an idea of the typical B-band radii and luminosities 
of normal galaxies at different redshifts. In order to compare with 
the models, in Fig. |2] the sample is split in three luminosity ranges 
indicated in the figure caption. The luminosities were calculated 
from the absolute magnitudes reported in the paper, and both Lb 
and Rb were rescaled to /i = 0.7 (they used h — 0.65). 

In Fig.[3l the sample is split in three z ranges indicated in the 
figure caption. As seen in Figs. |2] and [3] the few observed galaxies 
seem to be consistent, within the large scatter, with the radius evo- 
lution of our models. If any, the observed galaxies at a given z and 
radius, seem to be on average brighter in the B band than models, 
specially for higher redshifts. Since Lb traces recent SF rate, such 
a result may be related to the so-called downsizing problem. Such 
difficulty appears also in the similarity between the slopes of the 
Mt-R* and Lb-Rb relations showed by the models. From obser- 
vations, this slope is significantly steeper for the Lb-Rb relation 
than for the hh-R* one (e.g., Avila-Reese et al. 2008), reflecting 
this basically the fact that our models do not reproduce the observed 
correlation between color and mass or luminosity. While the pre- 
dicted slope in the B band seems to be only slightly shallower than 
observations, it is the slope of the R^-M^ relation that is steeper 
than the observational inferences presented in Avila-Reese et al. 
(2008; see also Pizagno et al. 2005). This difference is even more 
dramatic at lower masses when comparing models to the observa- 
tional results by Shen et al. (2003). In any case, we do not attempt 
to carry out a more detailed comparison of models and observations 
in the A/,-_R, ot Lb-Rb diagrams because these relations are actu- 
ally highly scattered and segregated by SB and morphological type 
(see e.g., Graham & Worley 2008). 

In order to infer reliable estimates of the average disc radius 
change with z, samples at different redshift bins, complete (cor- 
rected by the selection effects) down to a given luminosity, are nec- 
essary. The changes of the luminosity and mass-size relations in 
the range 1 S z < 3 of luminous (massive) field galaxy populations 
have been inferred recently from the Faint InfraRed Extragalac- 
tic Survey (FIRES) by Trujillo et al. (2006; see also Trujillo et al. 
2004). Due to the large number of objects in their sample, the au- 
thors are able to apply a completeness function to three redshift 
bins. Late-type galaxies are selected as those with the Sersic in- 
dex n S 2. Combining the analysis of FIRES data with the results 
obtained by using the Galaxy Evolution from Morphologies and 
SEDs (GEMS) survey at 0.2 < z < 1 (Barden et al. 2005) and 
tying both to the z ^ Q results from SDSS (Shen et al. 2003), Tru- 
jillo et al. (2006) presented a comprehensive picture of the change 
with z of the mean radius and its dispersion for galaxies more lu- 
minous than Lv ~ 3.4 W^^Lvq (rest-frame) or more massive 
than M« « 3 lO^^M© at any epoch. Their results for the n < 2 
(late-type) galaxies are reproduced in the upper (mass) and lower 
(luminosity) panels of Fig. |4j the error bars are the 2a scatter on 
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Figure 4. Upper panel: Crosses are the mean stellar radii normalized to the 
present epoch for observed disc galaxies more massive at each z than Af, 
3 1O^"/i^q'^M0 (from Tnjjillo et al. 2006). The error bars are the 2a scatter 
on the mean values estimated from the Log[_R*(z)/iJ»(SDSS)] distribution. 
The solid line is the /?* change normalized to z = of model galaxies 
that at each z have M, = 3 IO^^Mq (the behavior of the normalized 
Rt change for larger values of M,=constant. is quite similar to the one 
plotted here). Because the models correspond only to the "central" cases, 
there is not any scatter to report. The dotted line is the normalized halo 
radius change with z at a fixed A/h for all redshifts (eq.[T). Lower panel: 
Same as in the upper panel but in the V band, for galaxies more luminous 
at each z than Ly ~ 3.4 lO^^h^ffLvQ (rest-frame). The model galaxies 
are for Lv = 3 lO^^Lv-,© (solid line) and Lb = 2 W^°Lb,Q) (dashed 
line) at each z. 



the mean values estimated from the Log[i?e(z)/i?e('S'-DS'5')] dis- 
tribution. 

The solid lines in Fig.|4]show how the normalized stellar and V 
band effective disc radii change with z for models that at any z have 
M, = 3 lO^^Mo (upper panel) and rest-frame = 3 lO^^Lv,© 
(lower panel), respectively. In fact, the curves are quite similar for 
Af,=constant in the range IO^'Mq — 3 IO'^^Mq and Lv'=constant. 
in the range 2 10"'Lv,q - lO^L^,© (see Fig.^l. Thus, Fig.|4]cap- 
tures the typical shape of the normalized radius change with z for 
galaxies of the same mass or luminosity at different z's. Some ob- 
servations suggest that the normalized i?« change with z could be 
slightly steeper for the largest masses (Barden et al. 2005; Trujillo 
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et al. 2006; Franx et al. 2008). However, any possible dependence 
with mass is within the large scatter of the Rt-M, correlations at 
different z's and the observational uncertainties. Besides, at higher 
z's the samples are probably biased to observe earlier- type galaxies 
since they strongly emit in UV at their early formation phases; and 
early-type galaxies are more compact than late-type galaxies. 

The dotted line in the upper panel shows the redshfit change 
of the normalized halo virial radius, Rh given by eq. Our evo- 
lutionary models show that the relative radius change with z for 
a fixed disc stellar mass is much slower than in the case of dark 
matter haloes. For completeness, we also plot in the lower panel of 
Fig. |4l models in the B band, for rest-frame Lb ~ 2 10^^Lb,q 
(dashed line). The sequence is clear: the decrease of the normal- 
ized radius with z for a fixed luminosity (or stellar mass, which is 
closely related toa K band) is steeper as the pass-band is bluer. 

We did not include in our models the possible cooling flow 
regime for gas infall due to cooling times longer than the dynam- 
ical time in massive haloes. By truncating the late infall of high- 
angular momentum gas, the disc radius increasing is expected to 
be halted. Detailed calculations (e.g.. Button & van den Bosch 
2008) show that gas infall is truncated at z ~ for haloes of 
virial mass 14 ~ 6 lO^M© (corresponding to A/, f» lO^^M© 
and Lv ~ 7 lO^Lv'© in the case of the models presented here). 
As the halo is more massive than this value, the truncation of gas 
infall starts earlier. Therefore, the long cooling time in massive 
haloes is expected to work in the direction of flattening the late 
radius change with z of galaxies of masses A/. > 1O^°M0 (or 
Lv > 7 10^ Lv©), approaching even more models to observations 
in Fig. |4] We did not include either the possible feedback mass 
outflow mechanism, which introduces a decrease of /gai as Mh is 
smaller. Such a dependence is expected to produce a flattening of 
the Rt-M, relation at lower masses (Shen et al. 2003; Somerville 
et al. 2008), improving perhaps the comparison with observations 
in that regards this relation. 



4 CONCLUSIONS 

We have presented results on disc size change with redshift ob- 
tained by means of galaxy evolutionary models within the context 
of the ACDM cosmogony. Under the assumption of detailed angu- 
lar momentum conservation, discs are formed in centrifugal equi- 
librium inside growing ACDM virialized haloes, taking into ac- 
count the central halo (adiabatic) contraction due to the disc gravi- 
tational drag. The gas infall rate is assumed to be a constant fraction 
/gal of the halo mass accretion rate. The hydrodynamic equations 
for the gaseous and stellar discs are resolved coupled under the as- 
sumption of azimuthal symmetry and vertical virial equilibrium. 
The SF is triggered by the Toomre instability criterion and self- 
regulated by an energy balance between SN energy input and ISM 
turbulent dissipation. The driving parameters of the models are the 
halo MAH, the spin parameter \h, and the mass galaxy fraction, 

/gal- 

Since our main goal was to explore the general behavior of 
disc size evolution, here we took into account only "central" models 
constructed with the average MAH corresponding to a given Afh, 
a value of = 0.035 close to its mean, and a global average 
value for /gai=0.03; both \h and /gai are further assumed to be 
constant in time. We avoided the introduction of some (uncertain) 
astrophysical processes in order to explore in a transparent fashion 
the direct consequences of the ACDM cosmogony on disc size (and 
mass/luminosity) evolution; the assumed value of /gai is based on 



the success of the models to produce nearly flat rotation curves and 
realistic disc-to-halo velocity ratios. Our main conclusions are: 

- The model stellar and B luminosity disc effective radii sig- 
nificantly decrease with z (inside-out disc formation), the decreas- 
ing at high jz's being slightly steeper for the more massive galaxies 
(Fig.|2}. The radius growth in the B band is only slightly steeper 
than the one of the stellar disc, and at 2 = 0, Rb is 1.4-1.6 times 
larger than i?, . 

- The change of _R, with 2 at a fixed Af, -this is the kind of 
information able to infer from the observations- is much weaker 
than the radius evolution of the individual models (Fig. For 
z ^ 2.5, -R,(A/,=const)oc (1 + z)~^'^ for all the masses stud- 
ied here. This law is much shallower than the one corresponding to 
the haloes, which was widely used in the literature for comparisons 
with observations. On the other hand, in the B band, the decrease of 
Rb{Lb =const) with z is strong and it resembles the radius decreas- 
ing law of the individual models. For 2<0.75, i?_B(L_B=const)(x 
(1 + ^)-°■*^ andfor0.75<z^ 3, Rb(L B=const)(x (1 + 2)"^ ^ 
Since the radius growth of the stellar and B band discs is similar 
for the individual models, the difference in the radius change with 
2 for models of the same Af, and Ls at all 2's is basically due to 
the differences in the individual evolutionary paths of Af« and L b 

- The slopes of the _R,-A/, and Rb-Lb relations remain 
practically the same since 2 = up to 2 = 3; these slopes at 
2 = are 0.38 and 0.42, respectively (Fig. [3). Therefore, the 
changes with z of the relations zero-points, 7?* (2) / Af* (2)" '^** and 
Rb{z)/ Lb{z)'^'^'^ , describe basically the radius change curves at 
Af,=constant and Lb ^constant shown in Fig. [2] The similarity in 
the slopes of both relations disagrees with observational inferences 
at 2 = 0, which show a significant steepening from the stellar to 
the B band relation. Such a behavior is related to the well known 
color-magnitude (mass) relation, which our models do not repro- 
duce. 

- Model predictions are in reasonable agreement with observa- 
tional inferences of the typical R, and Rv change with 2 for disc 
galaxies more massive than 3 10"'Mq and more luminous than 
3.4 W^^LvQ (rest-frame) at any 2, respectively (Fig. |4ll. These 
large limits in A/« and Lv are dictated by the completeness func- 
tions required to introduce at different redshifts in order to elimi- 
nate the selection effects (Trujillo et al. 2006). In more detail, the 
predicted radius decrease with 2 for massive galaxies tends to be 
steeper than the one inferred from observations up to 2 ~ 1, spe- 
cially in the V band; for higher redshifts (up to the observational 
limit of 2 = 2.5 — 3), models and observations agree well. 

- Unfortunately, there are not yet available complete high- 
redshift galaxy samples that would allow to infer the typical stel- 
lar and luminous disc size evolution of sub-L« galaxies. Within 
the large scatters, our models seem to be consistent with small 
(not complete) samples of sub-i* disc galaxies with available rest- 
frame photometric properties (e.g., Tamm & Tenjes 2006) concern- 
ing the measured values of Rb and Lb at different 2's. 

The kind of model predictions presented here, unlike previ- 
ous ones made only on the level of galaxy populations at different 
redshifts, allow to understand the evolution of individual galaxies 
and its connection to the properties of the overall galaxy popula- 
tions. The model predictions are roughly consistent with the cur- 
rently available observational information on disc size evolution, 
showing that such an evolution is driven mainly by the ACDM halo 
evolution. The role of some astrophysical mechanisms in play (e.g, 
mass outflows, late gas infall truncation in massive haloes, stellar 
and AGN gas re-heating) seems to be of less relevance for the disc 
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size evolution. More observational inferences are required in order 
to constrain different models. 
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